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1. Introduction 



Thanks to recent algorithmic improvements and the availability of a new generation of ca- 
pability computers, simulations with dynamical Wilson fermions can now be extended to much 
hghter quark masses and bigger lattices. In recent years the QCDSF collaboration has significantly 
extended the number of ensembles of gauge configurations with Nf = 2 flavours of -improved 
Wilson fermions. These datasets cover a large range of quark masses as well as different lattice 
spacings and volumes. In principle, this puts us in the position to check for systematic errors that 
affect essentially all lattice calculations, i.e. finite size corrections, discretisation errors and errors 
due to extrapolations to the physical quark masses or the chiral limit. 

For the analysis of both, finite size effects and the quark mass dependence, results from chiral 
effective theories can be used to correct the simulation results or to guide extrapolations. Although 
the strategy of combining the results obtained from lattice simulations and chiral effective theories 
has led to a consistent picture for various observables, e.g. the pion decay constant [1], the mass of 
the nucleon [2], the axial coupling of the nucleon [3] or the nucleon electromagnetic form factors 
[4], this approach suffers from the paucity of lattice results in the region where chiral perturbation 
theory is expected to be applicable. 

Another source for systematic errors stems from the translation of the lattice results into phys- 
ical units. A quantity which is often used for this purpose is the Sommer scale tq. While this 
quantity has the advantage that on the lattice it can be determined with relatively high statistical 
accuracy, its experimental value is far less well known. To avoid depending on this experimental 
value, which relies on model assumptions, one can alternatively request lattice results for a quantity 
X, which has the dimension of a mass (e.g. the pion decay constant or the mass of the nucleon), to 
be equal to the con^esponding experimental result, i.e. Tq^'X'^' = roX'^'^P. The results from differ- 
ent collaborations suggest that the resulting value for ro is significantly smaller than the typically 
quoted experimental number ro ~ 0.5 fm. In this paper we will use ro = 0.467 fm which allows 
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for easy comparison with our previous results. Our most recent results for the nucleon mass (see 
next section and Fig. 1) and the pion decay constant (see [1]) suggest that the actual value is even 
smaller. To define our scale independent of the quark mass we extrapolate the measured values of 
^^'(jS, K"^) to the chiral limit using the ansatz 

lat,c 

In = Ao(j8) +A2(jS) (amps(j8, JC^))' , (1.1) 
a 

where A;(jS) =A;o+A,i(j8 - jSq) +A,-2(j8 - jSq)^ and j3o = 5.29. We fit our data in the range amps < 
0.5 and find x'^/^d.o.f. = 6.7/12 (see Fig. 1). 



2. Masses 



The mass of the nucleon, m^, and delta, wta, are experimentally well-determined quantities. 
However, for both quantities calculations using baiyon chiral perturbation theory (QxPT) suggest 
a rather non-trivial quark mass dependence. In an infinite volume wjn depends on the quai'k mass 
in the following way [2] : 



w^n("^ps) = Mq — 4ciwi 
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Even extending our fit interval to < mps < 650 MeV does not allow us to determine all parameters. 
We therefore restrict the set of free fit parameters to the nucleon mass in the chiral limit Mq, the not 
very well known low-energy constant (LEC) ci and the counter-term e[{X) (we use A = IGeV). 
For the other parameters, i.e. the LECs C2 and C3, the pion decay constant Fo and the nucleon axial 
coupling gAfi, we use the phenomenological estimates listed in Table 1. Our nucleon mass data and 
the fit ai^e shown in Fig. 2. Note that all our results seem to fall on a universal curve, indicating 
discretisation effects to be small. We therefore ignored 0{a-^) effects in our fit ansatz. We observe 
that mN('Mps = tuji) is consistent with experiment. Furthermore, we find c\ = —1.02(7) GeV^\ 
a value which is consistent with other estimates [11]. We have repeated our analysis for different 
values for tq. The results are shown in Fig. 1. If we set hjn = rn^^ we obtain ro = 0.457(3) fm. 

To extrapolate our results for the mass of the delta mA(mps) we fit these to the second order 
small scale expansion (SSE) expression [10] 
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Table 1: Phenomenological values used as input to our fits. 
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Figure 2: Lattice results for ton (left) and mA (right) together with a fit to Eq. (2.1) and Eq. (2.2). The open 
symbols indicate data which has not been included in the fit. The filled squares symbols in the left plot show 
preliminary data on larger lattices which has not been included in this analysis. 

where we fixed Fq to the value given in Table 1. We find that mA(?wps = mji) is equal to the ex- 
perimental value within statistical errors. The other fit parameters are consistent with expectations 
[10]: ai = -0.8(3) GeV"^ ~ci and /ia = 1.5(4) < 9^a/5. 

3. Nucleon axial coupling 

We will now consider the form factor of the nucleon axial current Ga{Q^) at zero momentum.^ 
The axial coupling constant gA = Ga(0) is determined from the renormalised axial vector current 
= Za (1 +bAamq)A^, where aniq = {1/k—I/ k^) /2. Here we only consider the iso-vector case 
where contributions from disconnected terms cancel. While Za is known non-perturbatively [3], 
bA is only known perturbatively and is computed using tadpole improved one-loop perturbation 
theory. Note that for forward matrix elements there is no need for improvement of A^. 

The quark mass dependence of the iso-vector nucleon axial coupling gA{mps) has been calcu- 
lated using Heavy Baryon ;^PT (HB;^PT) [3]. These calculations have been performed both in the 
infinite volume limit as well as for a finite spatial cubic box of length L. We define gA (wips , L) = 
§a('?^ps) + AgA ('Wps , Z.) , where AgA{mps,L) denotes the finite size effects. 

We fit our data to the results obtained in [3] fixing the values for Fq, the leading axial NA 
coupling ca and the NA mass splitting at the physical point Aq to the phenomenological values 
given in Table 1. We furthermore set B2o('^)'^'^^ = and use A = m„. A comparison of our data 
and the resulting fit is shown in Fig. 3. 

The raw lattice data is significantly smaller than the experimental value. The fit to the HB;^PT 
expression suggests that this might be due to significant finite size effects and a rather strong quark 
mass dependence in the range < mps < 300 MeV. gA{nips = nij^) is found to be smaller than the 

'Results for non-zero momentum have been presented by W. Schroers at this conference [4]. 
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Figure 3: Lattice results for gA together with a fit to an expression obtained from HB;)fPT. The filled symbols 
show the lattice results after finite size effects have been corrected. 

experimental value, which is however not significant given the large statistical errors. The other fit 
parameters are of natural size. 

4. Moments of unpolarised structure functions 

Finally, we consider the lowest moment of the unpolarised nucleon structure functions, (x) = 
AI q{0), where q is the first moment of the PDF H'^{x, ^ , Q^) at = 0.^ This moment is deter- 
mined from the matrix element 



For O(fl') -improvement of the operator we use the perturbative results obtained in [14], where we 
inserted the boosted coupling constant (g*)^ = E^/uq. The unknown improvement coefficients are 
set to zero as the corresponding improvement operators turn out to be small. The renormalisation 
is both scale and scheme dependent. We have calculated the renormalisation constant AZj,2b(fl'), 
which translates our lattice results into RGI, non-perturbatively using the RF-MOM method. For 
comparison with other results we have to convert our numbers to MS. The corresponding factor 



The lattice results for {x) in the iso-vector case have been found to be significantly larger than 
the experimental value. It had been suggested that this quantity may become much smaller at very 
light quark masses [15]. This has been confirmed by recent calculations in the framework of baryon 
chiral perturbation theory (B;^PT) [8] . 

We compare our lattice results with these calculations both in the iso-vector case (u — d) and 
the iso-scalar case {u + d). Note that the latter might be affected by contributions from disconnected 
terms which we have not calculated so far. In our fits we use phenomenological input for the pion 
decay constant Fq, the nucleon mass Mq and the nucleon axial coupUng g^.o as listed in Table 1. 

^Results for A| q at 2^ > have been presented by M. Ohtani at this conference [12]. 



{N{p)\ uy^^'D^'K \N{p))=2Alo[p^"p^']. 



(4.1) 



AZ^^(At = 2GeV) 
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is calculated perturbatively using A^^ = 261(17)(26) MeV [13]. 
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Figure 4: The left plot shows our results for (x)("^^)(mps) at j3 = 5.29 (the results for other values of 
j3 are similar). The solid line denotes a fit to an ansatz obtained from BjPT [8]. The right plot shows 
(x) (»zps = OT;t) as a function of the lattice spacing together with a fit to a constant. 

As suggested in [8] we assume the coupling Aajo — (Ax)'"^'^'). We finally end up with two free 
fit parameters, i.e. a" and c^{X = \ GeV) as well as {d^^) and c<){X = \ GeV) in the iso-vector 
and iso-scalai- case, respectively. With the statistical eiTors and the number of fit parameters being 
sufficiently small we fit the data for different values of j8 separately restricting the fit range to 
< nips < 650 MeV (see Fig. 4 and 5). 

In the iso-vector case we still find little indication for (jc)("^'^)(?Mps) becoming smaller for 
mps — > mji. Currently this is not in contradiction to the ^X?T results since quark masses may 
still be too large. Furthermore, results at light quark masses may be affected by finite size effects. 
Within statistical eiTors the discretisation effects seem to be small. In the iso-scalar case we find 
evidence for (;c)("+'^)(mps) to become smaller at lighter quark masses. Our results for nips = m^i 
are close to the phenomenological value. 



5. Conclusions 

In this paper we presented recent progress in improving our control on the systematic errors 
for various quantities in the baryon sector. Finite size effects were found to be small in the case 
of the nucleon mass win and large in the case of the nucleon axial coupling gA- Investigation of 
such effects for the delta mass mA and the lowest moment of the unpolarised nucleon structure 
functions {x) is still pending. For all quantities considered here discretisation effects seem to be 
small, although it remains difficult to distinguish potential 0{a^) effects from other systematic 
eiTors. The largest uncertainties stem from the extrapolation of our results to the physical quark 
masses. We utilised results obtained in xPT to perform these extrapolations. Although most of our 
lattice simulations have been performed at quark masses outside the region where ;^PT is expected 
to converge, our fit parameters turned out to be of natural size and consistent with what is expected 
from phenomenology. 

To get better control on the quark mass dependence, simulations at significantly lighter quark 
mass are mandatory. Fortunately, with recent algorithmic improvements and even faster computers 
becoming available exploring the quark mass region with mps < 300 MeV starts to become feasible. 
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Figure 5: Same as Fig. 4 but for the iso-scalar case 
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